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Abstract 

We present a method for calculating the photoproduction of jets at HERA 
based on Next to Leading Logarithm QCD calculations. It is implemented 
in a Monte-Carlo generator which allows us to easily compute any infra-red 
safe cross sections for 1 or 2 jet observables using various jet reconstruction 
algorithms. We focus on the possibility of extracting the gluon contents of 
the photon and of the proton from present and future HI and ZEUS data. 
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1 Introduction 



Electron-proton scattering at HERA is dominated by the exchange of a quasi-real 
photon and a fraction of j-p collisions leads to the production of high p T jets. 
Therefore these processes can be predicted by perturbative QCD and they can be 
used to investigate the structure of the photon in a way complementary to the study 
of the deep inelastic scattering e-j. In the latter case, one can measure directly 
the quark density inside the photon whereas the gluon density is constrained by the 
evolution equation. On the contrary, in photoproduction, one can probe directly the 
gluon contents of the photon. 

Leading Logarithm (LL) results have been available for a long time, but they 
are plagued by sizeable uncertainties coming from the dependence on unphysical 
scales. Therefore, Next to Leading Logarithm (NLL) calculations are essential to 
describe the photoproduction of jets. To reach this goal, we have adapted a method 
developped previously to deal with the production of two high-p^ hadrons in hadron 
collisions [1]. We apply it to build a Monte-Carlo generator which is able to produce 
a set of partonic events on which one can apply any jet reconstruction algorithm to 
produce a set of jet events. With the latter we can easily compute any infrared safe 
cross sections. This technique gives a lot of flexibility to the study of various jet 
algorithms and observables. 

In particular we will apply it to the study of two jet observables in which the 
distribution function variables are well constrained by kinematics. This allows us to 
disentangle the distribution functions from the hard subprocess more easily, because 
less convolutions are involved in the calculation of the cross sections. Some results 
obtained with this approach have already been reported in ref. [0. 

The paper is organised as follows. The theoretical framework is described in 
section 0. We also compare our approach with those of other authors. In section |3| 
we present some applications of our work to 1-jet cross sections and we compare 
our results with those already obtained with an analytical approach. Section |] is 
devoted to the study of 2-jet observables. We will examine whether it is possible to 
accurately measure the distribution functions of the gluon in the photon and in the 
proton. In section [5] we study recent HI and ZEUS data and the constraints they 
put on the gluon distributions. The conclusions are in section |. 

2 Description of the method 

To calculate cross sections and isolate collinear and infrared singularities, we used 
a "phase space slicing method" . We start from the 2 — > 3 partonic squared matrix 
elements and virtual corrections evaluated by Ellis and Sexton J3|. Collinear and 
infra-red singularities lead to poles in 1/e and 1/e 2 when using dimensional regular- 
isation. For a generic real process 1 + 2^3 + 4 + 5, at least two partons have high 
transverse energy E? (3 and 4) and only one can be soft (5). In order to extract 
the singularities, we cut the phase space in several parts: part I where E?5 is less 
than a given scale prm and part II where Et$ > PTm- Part II is divided in three 
parts: Ha (resp. lib) where parton 5 is within a cone around parton 3 (resp. 4) 
called C3 (resp. C4), lie where parton 5 is outside C3 and C4. Parton 5 is in Ci 



if ((05 — (pi) 2 + (775 — rji) 2 )^ < R c . Here r\ = — logtan(#/2) is the pseudo-rapidity 
and is the azimuthal angle. Part I contains infra-red singularities and collinear 
singularities in the initial state and parts Ha and lib contain collinear singularities 
in the final state whereas part lie is finite. 

The contributions of regions I, Ha and lib are calculated analytically and the 
infra-red singularities are cancelled by the corresponding ones in the virtual terms. 
Initial collinear singularities are factorised in the parton distributions and the final 
collinear singularities disappear when integrating on the relative momentum between 
parton 5 and the parton with which it is collinear due to energy momentum sum 
rules. 

The finite parts remaining after the cancellation of singularities have been ana- 
lytically computed using Maple [f|. Large logarithms logpr m , log 2 pTm and logi? c 
appear in the collinear and infrared contributions. They are cancelled by similar 
terms from part lie so that the total cross section is independent of these unphys- 
ical cuts. It should be noted that we have kept only the logarithmic terms in the 
calculation of contributions I, Ha and lib, neglecting terms of order 0{pTm logprm), 
0(R 2 log PTm), and less singular terms. 

Using the Monte-Carlo package BASES ||, our program generates quasi 2^2 
events corresponding to collinear, Born and virtual contributions (we take pTm <C 
E™™ and R c -C 1) and 2 — * 3 events corresponding to part lie. For the latter, a 
jet reconstruction algorithm is applied. Finally, these events are histogrammed || 
in order to give any cross sections we are interested in. We have checked that 
the cross sections do not depend on the unphysical cuts pTm and R c in the region 
where 0.005 < pT m < 0.1 and 0.01 < R c < 0.1. From a numerical point of view, 
the compensation between this cut dependence arises mainly between positive real 
contribution of part lie and negative contributions of part I, Ha, lib and virtual 
corrections ; we need to generate a sufficient number of events in order to obtain 
a small error after the compensation. We take the greatest possible values for the 
cuts PTm and R c (i.e. pTm = 0.1 GeV and R c = 0.1) in order to lower the size of 
this compensation which is then typically of the order of 1 for 5. 

This approach is applied to the direct and resolved (proportional to the photon 
distribution functions) parts of the cross section. However these contributions de- 
pend on the convention adopted in the subtraction of the collinear singularities. Un- 
less explicitly specified, we use the MS factorization and renormalization schemes. 

Several authors have used similar approaches in the calculation of jet-photopro- 
duction cross sections. Harris and Owens also built an event generator ; but their 
phase slicing is based on cuts put on E5 (in the parton CM) and on the Mandelstam 
variables of the subprocess which control the collinear divergences. Klasen and 
Kramer || introduced a single invariant cut on the Mandelstam variables. On the 
other hand Frixione and Ridolfi [[| made their calculations with the subtraction 
method ||10||. We made several numerical comparisons with the results of Klasen 



and Kramer, and found a good agreement. 

Finally let us note that it is very convenient to use variables Et, T] and <fi to 
define regions of the phase space in which we perform analytical calculations. Indeed 
these are the variables commonly used by experimentalists to define jets and cuts 
in their phase space. If pT m and R c are smaller than respectively the lower bound 



of the transverse momenta and the width of jets, we never integrate (thus defining 
an inclusive measurement) in a phase space region in which the experimentalists 
perform an exclusive measurement. Since these experimental bounds are much larger 
than the domain in which our neglecting non logarithmic terms is valid this is in 
fact not a constraint at all. On the contrary with the invariant cut method, this can 
only be achieved with very small cuts. 



3 Single jets 



In this section we present some numerical results for one-jet cross sections in order to 
compare our predictions with those obtained in a preceding paper pTTj] . Our inputs 
are the following. At the HERA collider (y^ = 300 GeV), electrons produce photons 
with small virtuality Q 2 . We use the kinematical conditions of the ZEUS collabora- 
0: 



tion 



Q 



2 

max 



1 GeV 2 and 0.2 < y < 0.85 where y = E^/E e . The spectrum of 



the quasi-real photon is approximated by the Weizsacker- Williams formula. 



F i( v ) = ^S I + O^zjH log QL*(i-v) _ 2(1 -y) 

2tt \ y * m 2 y 2 y 

For the proton distributions we take the CTEQ4M parametrization [O and 



for the photon distributions the GRV parametrization [14| transformed to the MS 



scheme. We use five flavours, = 296 MeV, a renormalization scale a and a 

' MS ' ^ 

factorization scale M equal to the transverse energy Et of the jet. For all our 
calculations, we use for a s (/i) an exact solution of the two-loop renormalization 
group equation, and not an expansion in log^-. The cross sections are higher by 
some 2.5 % when the exact a s (fi) is used with ji ~ 15 GeV. 

Here we present results only for the resolved contribution. Jets are defined with 
the fc r -algorithm [HJ. We compare our results with those of a previous analytical 



calculation [[llj, for the transverse momentum distribution in Fig. |l]a and for the 
rapidity distribution in Fig. 0b. We can see that the agreement is quite good between 
these two sets of results. Similar comparisons hold for the direct part. 

We do not pursue the study of 1-jet cross sections, because they do not constrain 
the parton distributions as well as the 2-jet cross sections do ; the latter offer more 
kinematical possibilities of control of the x-variable of these distributions. 



4 Dijets cross sections and the photon structure 
function 

In this section we study the dijet cross sections and put the emphasis on variables 
and cross sections which give access to the gluon distribution functions. Comparison 
between data and theory is postponed to section 5. The dijet cross section, as a 
function of the transverse energy Et3 and the jet rapidities 7/3 and 774, is given by a 
product if the subprocess is a 2 — > 2 reaction (LL approximation) 
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Figure 1: a) Single jet resolved cross section da/dE T integrated in the rapidity range 
.5 < 7] < 1.5 (histogram) and compared to the analytical result of b) da/dt] for 
the resolved contribution with 20 GeV < Ex < 21 GeV. 



(h - M) x p F p 6 (x p) M) ^" (/i) , (2) 



dEindi] 3 dr]4 ^ e e v e ' ' ' ' dt 

± ° a,b,c,d 

da/dt is the a + 6 — > c + d cross section, and F£(Fp) the parton distributions in the 
electron (proton) (ac p = E T3 (e m + e r > 4 )/2E p and x e = E T3 (e~ m + e~^*)/2E e ) ; M 
is the factorization scale and \i the renormalization scale. The direct contribution 
corresponds to a = 7 and _F e 7 is the Weizsacker- Williams formula ([J). In the resolved 
case, _F e a is given by a convolution of _F e 7 with the parton distribution in the photon 



F?{x e , M) — [ dy dx 7 F2{y) F^, M) 5{x lV - x e 
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(3) 



If the photon energy E 1 = yE e is known (for instance by tagging the outgoing 
electron), we can measure x 7 and define the observable 



x 7 F 7 1 (x 7 ,M) / dE^ 3 dr) 3 dr]4 5 ( x 
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Thus the dijet cross section written as a function of x 7 is proportional to the parton 
distribution in the photon F"(x 7 ). One observes that the direct contribution, with 
F" = <5 7a 5(l — x 7 ), leads to a peak in the cross section at x 1 — 1. 
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Figure 2: A comparison between the distributions da/dx^ L and da/dx° bs , with the 
following cuts on the two high E T jets: E l ° adin9 > 12 GeV and E^ ailin9 > 10 GeV 
and < T) < 1, and the fc^-algorithm for the definition of jets. 



NLL QCD corrections to the LL expression (Q) blur its simple kinematics. There 
are contributions with three jets in the final state, and we can no longer fix x 1 = 



Et 2e Vl > because the third jet is not observed. However we can follow the strategy 

i=3 

of the ZEUS collaboration which defines the variable (jets 3 and 4 are the jets 
with the highest transverse energies) 
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(E T3 e~ m + E Ti e 
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(5) 



and observe the dijet cross section da/dx° bs . However this definition of x° bs may 
lead to infrared sensitive cross sections. Indeed fixing Et3 and E Ti strongly con- 
strains the available phase space of the unobserved partons (for instance parton 5 
in section 2) and forbids a complete compensation between real and virtual NLL 
corrections. This results in cross sections containing, for instance, terms propor- 
tional to log (l - (E™ n /E T3 ) 2 ) after an integration over Eta from a lower bound 

E™™ (smaller than Etz) has been performed. The cross section is not defined at 
Et3 = E™™, although it is integrable. Therefore the variable x° bs defined in (|5D 
should only be used with cross sections integrated over E T 3 and E T 4 ; moreover the 
integration range of E T 3 and E T i should not have the same bounds. A discussion of 
this condition may be found in ref. 0. 
For these reasons, we prefer the variable 
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LL 
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(6) 



which depends on the transverse energy of only one jet. It can be associated with 



cross sections in which the energy of the second jet is not measured. One observes 
that x^ L may take values larger than 1.0. 

One must also observe that the cross sections da/dx^ L or dajdx obs are singular 
when x^ L or x obs approaches 1. Indeed we have the condition 

E T3 e-^ + E Ti + E T5 
(the sign = corresponds to the direct case), which means 



1 - x obs > 



T5 e 



2E,, 



(7) 



When x° goes to 1, the phase space of parton 5 is severely restricted. This results 
in log(l — x° bs ) terms generated by the NLL corrections. We obtain a similar result 
with x^ L —> 1, although the phase space of parton 5 is less severely constrained than 
by condition (0). Therefore we expect a smoother behaviour of da/dx^ L . This point 
can be verified in Fig. ^| in which we display da jdx obs and da/dx^ L for the direct 
term. Besides the region very close to 1.0, the two distributions are very similar Q. 
(The AFG [|TJ and ABFOW @ distributions with N f = 4 have been used for this 
calculation). Currently experimentalists integrate x° bs over the range .75 < x obs < 1. 
Therefore the singular behaviour of the cross section is smoothed and should not 
forbid the phenonemological application of NLL calculations. 

We can also see from Fig. |2] how much NLL corrections modify the x 7 distribution 
which is proportional to 5(1 — x 7 ) at the LL accuracy. NLL corrections generate 
terms with x^ s or xz different from 1 and the simple picture of a photon directly 

interacting with a quark of the hard subprocess at OC'y X IS lost. But one must 

also keep in mind that the separation of the cross section da jdx obs into a direct part 
and a resolved part is factorization scheme dependent and that da dvrect / dx obs has no 
physical sense on its own (the same remark is valid for x^ L ). 

To study this problem we computed the observable da/dx^ L in two factorization 
schemes, the so-called DIS^ and MS ones. For the resolved part, we used for the 
former case the GRV distributions of quarks inside the photon and for the latter the 
AFG distributions. (Here we also use the ABFOW proton distributions). We can 
see that both the direct and resolved part are very different in these two schemes 
(Fig. |^), but that this difference is much smaller for the total cross section which 
is factorization scheme independent. The remaining difference partly comes from 
different hadronic inputs in the two sets of distributions. 

We now turn to the study of the gluon contents of the photon and to the pos- 
sibility of constraining it through the observation of the distribution da/dx obs . The 
sensitivity of da jdx obs to the gluon distribution does not only depend on the value 
of x° bs , but also on the various kinematical cuts imposed on the 2-jet phase space. 

4 Note that the shape of the cross sections around = 1 depends on the width of the x 7 -bin, 
because the Born and virtual contributions are proportional to 6(1 — a; 7 ). If the width is too small, 
the cross section may even be negative (the positive real contributions do not compensate anymore 
the negative virtual contributions). 
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Figure 3: Distribution da/dx^ L for the direct, resolved and total contributions with 
the AFG (MS scheme) and GRV (DIS 7 scheme) photon density as well as ABFOW 
proton density, with the following cuts on the two high E T jets: E T > 9.2 GeV and 
— 1 < i] < 2, and using a cone algorithm with R = 1 and no R sep to define the jets. 
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with the following cuts for the two high E T jets: 



Figure 4: Cross sections da/dx° bs 

Reading > ^ Qe y ^ ^railing > 1Q Qe y ^ 2 different cutg Qn their rapidit i eS) 

and using the /^-algorithm to define the jets. Curves with the gluon distribution 
varied by ±30% around the AFG one are also displayed. 



Clearly the region of positive and large rapidities corresponds to small values of 
x obs and to an enhancement of da/dx° bs . This is verified in fig. 4. We divided the 
rapidity range in the lab frame for the two leading jets in subintervals: < rj < 1 
and 1 < rj < 2. We have also used asymmetric cuts on the transverse momenta, as 
it is done by experimental collaborations, in order to avoid the instabilities which 
appear when these momenta become close to each other. Then in each of these 
rapidity intervals we compared the x° bs distribution obtained with the AFG photon 
density with cross-sections for which we have artificially reduced and increased the 
gluon distribution by 30 %. We found that the influence of the gluon increases with 
the rapidities and in fact the differences between these three curves become sizeable 
only when 1 < rj < 2. Indeed an increase of 30 % of the gluon density results in an 
increase of approximately 25 % for the cross-section around x° bs = 0.2. Therefore a 
determination of the gluon contents of the photon would require to use such cuts on 
rapidities. Such a study would be able to test the gluon density Fj>(x) in the region 
0.2 where it is very poorly known. 



5 Comparison with HI and ZEUS data 

In this section we analyse recent HI and ZEUS data to assess the possibility to 
put constraints on the gluon distributions in the photon and in the proton. With 
this intention we investigate the sensitivity of various cross sections to changes in 
the gluon distributions, and we compare their variations to the experimental errors. 
In this way we obtained an estimate of the accuracy with which gluon distributions 
can be extracted from present data, and from future high statistics experiments. We 
modify the gluon distributions by increasing or decreasing their normalizations by a 
few tens of percents. This method has the advantage of leaving unchanged the well- 
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Figure 5: The cross section da/drj2 compared with ZEUS data [P^j . (Full error 
bars : statistical errors ; dashed error bars : systematic errors including energy- 
calibration errors). Full line: .5 < y < .85. Dotted lines: .45 < y < .85 (upper) 
and .55 < y < .85 (lower). The dashed line is obtained with the distribution g/7 
increased by 20 %. 



determined quark distributions and to quantify the gluon modifications in a simple 
way. Of course gluon distributions are constrained by other experiments. But, as 
discussed in the introduction, the gluon in the photon is not well constrained by DIS 
77* data. As for the gluon in the proton we are going to study observables sensitive 
to the distribution at small x ~ 0.02. In this region the gluon is pinpointed by the 
slope of F 2 and it was shown by the CTEQ collaboration that a variation of about 
±10% of Fg in this region and at a scale of 100 GeV would cause clear disagreements 
with DIS plus Drell-Yan data ]2B|]. However this conclusion was not obtained by 
performing an error analysis but only by tuning the gluon parametrisation, which 
might artificially reduce the range of variation by being too restrictive. A more 
recent study obtained a gluon density 30% bigger than CTEQ4M about x ~ 0.02 
and at a scale of 20 GeV [ 27| . Errors on the gluon determination are also presented 



in ref. [01. Thus we think that it is interesting to find how much photoproduction 



can constrain this density in this x-range. 

The interested reader may find more global comparisons between theory and data 
in ref p0| , ^2|, j23| . The inputs we use in this section have been defined in section 



3 where we computed single-jet cross sections. More details on the kinematical 
parameters used for dijet cross sections are given belowQ. 

5 To follow HI and ZEUS conventions, we call jet 1 and jet 2 the jets with the highest energies ; 
we used the labels 3 and 4 in the preceeding sections. 



The ZEUS collaboration presents cross sections which depend on the jet trans- 
verse momenta, or on the jet rapidities p3 |. Therefore we do not have a direct access 



to da/dx° bs , however we can study data corresponding to regions of the jet phase 
space in which the role of the gluon in the photon or that of the gluon in the proton 
is enhanced. These regions are essentially defined by large or small values of the jet 
rapidities: large rapidities correspond to small x 1 and small rapidities to small x p 
(see the LL expressions for x p and x e below formula @))- 

Let us start with the case of large jet rapidities. Experiment and theory are 
compared in Fig. |5|. The dijet cross sections in this figure correspond to events with 
at least one jet with transverse energy larger than 14 GeV, the transverse energy 
of the other jet being larger than 11 GeV ; it is integrated over 1 < rji < 2. The 
real photon kinematical domain is specified by Q 2 max = 1 GeV 2 and .5 < y < .85 (cf 
expression (P). The jets are defined with the /^-algorithm |H| and the scales M 
and fj, are set equal to the transverse energy of the most energetic jet. 

To test the sensitivity of da/drj2 to the gluon density in the photon, we increased 
the gluon distribution uniformly by a factor of 1.2. Comparing the full and dashed 
curves, we see that this factor produces an increase of the cross section for large 
values of 772 by some 10 %. As expected the backward region, corresponding to large 
values of x 7 , is not affected ; the gluon distribution in the photon decreases much 
faster than the quark distributions. In this region da/dr/2 rapidly varies with r] 2 and 
better comparison with data is obtained by integrating da/dr]2 in the experimental 
bins. This has been done for the two bins: —1. < 7/2 < —.5 and —.5 < r/2 < .0. The 
agreement between theory and data is quite good in the second bin but the data at 
larger 7/2 favor a larger gluon distribution in the proton. But one must notice that 
the present experimental errors are larger than the effect produced by an increase of 
the gluon distribution by 20 %. In the backward region, a clear discrepancy appears 
which could be attributed to hadronization effects as discussed in ref. ||24|| . 

It is also worth again noting the great sensitivity of the cross section to the y- 
range of the photon. The energy of the incoming photon is reconstructed from the 
final hadron energies with the Jacquet-Blondel method. Various corrections have to 
be applied to this "photon energy" y JB in order to obtain the true photon energy y 
2"5| . In Fig. ^ we show the effect of a 10 % error on the determination of the lower 



bound of the y variable. It is very large for negative values of 7/2. Even in the forward 
region, this 10 % error produces an effect much larger than the 20 % variation of 
the gluon distribution. Actually this "error" is included in the discussion of the 
systematic errors quoted in p5| . But because the dispersion of y JB around y may 
reach 10 %, a better theoretical prediction could be obtained by taking into account 
the dispersion of the upper and lower bounds of the variable y. 

In Figure |^ we present the results of a similar study done for negative values of 
rji. In this kinematical range, we can infer from fig. [|that the cross section depends 
very little on moderate changes of the gluon distribution in the photon. Therefore 
this kinematical range is well-suited for a study of the gluon in the proton at x p ~ 
2E T /2E p ~ 15/880 = .018. Here we choose data corresponding to the following 
?/-range: .2 < y < .85. Also displayed are predictions obtained for .2 < y < .8, and 
those obtained with a gluon distribution in the proton multiplied by a factor .8. We 
see that a small modification of the y-range produces an effect similar to a change 
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Figure 6: The cross section da jr\i compared with ZEUS data |23|]. (Full error bars : 
statistical errors ; dashed error bars : systematic errors including energy-calibration 
errors). Full line: .2 < y < .85. Dotted line: .2 < y < .80. The dashed line is 
obtained by reducing the distribution g/P by 20 %. 

of the gluon distribution by 20 %. Here again, the experimental errors are larger 
than the effects due to a modification of the theoretical inputs, but there is a slight 
indication that data prefer a smaller gluon distribution in the proton. 

Finally we turn to the -&r-spectrum obtained by the ZEUS collaboration. Here 
again we choose to study the large rapidity region with the dijet cross section inte- 
grated over the ranges 1 < 771 < 2 and 1 < r] 2 < 2: 



da/dE l ^ admg = J drji J di] 2 



da 



7 T-Xeadinq j -■ 

dE T M dr/i dr/2 



fading - g transverse energy of the leading jet (highest Et)- The transverse 
energy of the other jet is constrained by the condition 11 GeV < Et < E^ adm9 . Our 
results are shown in Figure 7. Predictions for the range .25 < y < .85 and for a gluon 
in the photon increased by 20 % are also shown. A slightly better agreement with 
data is obtained in the latter case (hardly distinguishable on a logarithmic plot). 

The HI collaboration presented results |20| under a form which is very close to 
the one advocated in the preceding section ; the authors chose to give the cross 
section da/{dx°^ s d Log((£^ s ) 2 /GeV 2 )) with the definitions 



4* = ET1 ± ET2 (9) 

where Eti and E?2 are the transverse energies of the two jets with the highest 
transverse energies in an event. Other kinematical requirements are 



100 



> 

O 

Oh 10 



<Q f-H 



t 1 r 



10 



1 <r| <2 



1 < T| < 2 

2 




.20<y<.85 



20 



30 



40 



50 



60 



E leading 
T (GeV) 



Figure 7: The cross section da / dE l ^ adm9 compared to ZEUS data f23[. (In this 
figure the energy-calibration errors are not included in the error bars). Full line: 
.20 < y < .85. Dotted line: .25 < y < .85. The dashed line is obtained by increasing 
the gluon distribution g/j by 20 %. 
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The cuts on the "real" photon variables are Q 2 max = 4 GeV 2 and .2 < y < .83. We use 
a cone algorithm |2l| with R = .7 and no R sep in agreement with ref. f20 |. Moreover 



we avoid double counting of jet configurations by choosing the jet of highest Et 
(made of two partons) when the parton configuration also allows to construct two 



jets (made of one parton each). The scales M and fi are set equal to E^ ts . 

In Fig. H we compare our results with HI data in the range 2.30 < Log I ) < 

2.50 where E^ ts is large enough to allow us to neglect (in a first approximation) 
hadronisation effects and underlying event contributions. A good agreement^ is 



Unlike other observables, here we disagree with the predictions of ref. []22|. Our values are 



o 



■s 




Figure 8: The cross section da/(dx° bs dLog((E^ it /GeV) 2 )) in the range 2.3 < 



Log(^ e VGeV) 2 < 2.5 PJ. Full curve : .2 < y < .83 ; dashed curve : .25 < y < .83. 



obtained with data, except for x° bs close to one (the "Direct" domain) where theory 
overshoots data. It must be again noted that the theoretical curve is very sensitive 
to the photon energy range y = E 7 /E e . A change of the lower limit from y = .20 to 
y = .25 leads to the dashed line in fig. |8| ; the cross section decreases for x° bs ~ 1.0. 
However it is unlikely that the disagreement for x° b ~ 1 between theory and HI 
data could be explained by a dispersion of the upper and low bounds in y. (The 
variation of y\ ow by 25 % studied above is certainly bigger than the experimental 
systematic error). A decrease of the gluon distribution in the proton would decrease 
da/dx° bs in all bins in x obs . It is compatible with data at low x° bs , but its effect 
at large x° bs is not sufficient to put theory in agreement with data. The resolved 
contribution is also important in the highest x° bs -bin, a kinematical region which 
explores the quark contents of the photon at large x 7 . In this domain the quark 
distribution is given by the pointlike component (non pointlike contributions of the 
"Vector Meson Dominance" type are negligible) which cannot be modified in an 
arbitrary way, and adjusted to data. For instance the difference between the AFG 
and GRV parametrization is small in this region (see fig. 3). 

The first conclusion that we can draw from this study of HI and ZEUS data is 
that overall there is a good agreement between experiment and theory. However the 
systematic errors are non-negligible and correspond roughly to variations of the the- 
oretical predictions coming from modifications by ± 20 % of the gluon distribution 
normalizations. Therefore it appears difficult to constrain the gluon distributions 
with an accuracy better than a few tens of per cents with the present data. 

Part of the errors should cancel in the ratio 

higher at small x° bs , by a few tens of percents ; this cannot be explained by the different values of 
Nf used in the calculations (4 flavors in ref. and 5 in this paper). 



which could allow a better determination of the x shape of the distributions func- 
tions. However we would thus loose any information on the absolute normalization. 



6 Conclusions 

In this paper we described a new NLL event generator for photoproduction reactions 
involving the direct and resolved contributions, and we used it to assess the possi- 
bility to measure the quark and gluon contents of the photon from photoproduction 
experiments. (Actually we assumed that the quark distributions in the photon were 
fixed from 77* DIS experiments and we concentrated on the gluon contents). 

The cross section da/dx^, where x 7 is related to the scaled momentum x 7 of par- 
tons in the photon, is quite appropriate to a study of these contents. We discussed 
two definitions of 00 ' /y . tC ■y OC ° , the well-known definition of the ZEUS collabora- 
tion, and variable which reduces problems of infrared sensitivity. Then 
we showed that the cross section da/dx° bs is sensitive to the gluon density of the 
photon only if we require the two jets with the highest transverse energy to have 
positive rapidities. 

A relatively good agreement between theory and experiment is found when con- 
fronting the predictions, obtained with the CTEQ4M and GRV distributions, with 
HI and ZEUS data. It is interesting to note that HI and ZEUS data are compatible 
with a 20 % increase of the gluon contents of the photon (Fig. 6, 771 ~ 1.5, Fig. 8, 
pleading ^ Q e v) and with a 20 % decrease of the gluon contents of the proton 
(Fig. 5, x° bs ~ 1 where the direct contribution is important ; Fig. 7). Unfortunately 
this remark cannot be made more quantitative, because the systematic errors, es- 
sentially coming from the uncertainties in the measurement of the jet energies, are 
large and of the same order as the variations of the theoretical predictions due to 
gluon distribution modifications. One must also keep in mind the great sensitivity 
of the cross section to the photon energy-range. The knowledge of the resolution 
of the variable y (obtained from the Jacquet-Blondel variable Ujb) should allow a 
more accurate prediction. However it is interesting to remark that 1996 and 1997 
preliminary ZEUS data also favor a larger gluon distribution in the photon. 

Until now we have not discussed the theoretical uncertainties coming from the 
scale dependence of the cross sections. These uncertainties were carefully studied 



in ref. || and in ref. |22| , and found to be of the order of a few tens of percents 
for large variations of the values of the factorization and renormalization scales. In 
this paper we continue this study by looking at the sensitivity of the cross section 
do I (dx" bs dLog(E^ ts ) 2 ) which can be compared with HI data (Fig. 5) ; we make this 
study for the range .75 < x° bs < 1 and with the kinematical conditions of Fig. 5. 
The scales are M = /j = kE^ s and k is varied between .5 and 2. Our results are 
summarized in table 1. 



K 


Direct contribution 


Resolved contribution 


Total contribution 


.5 


.634 


.384 


1.018 


.75 


.578 


.432 


1.010 


1.00 


.547 


.453 


1.000 


1.50 


.512 


.474 


.986 


2.00 


.491 


.482 


.973 



Table 1 : The sensitivity of the direct and resolved photoproduction cross-sections 
to the renormalization and factorization scale M = \x = k ■ E^ ts (normalized to the 
total contribution at k — 1.0). 



We note that the cross section da/dx° bs is quite stable and varies by less than 
5 % when the scales M 2 = /i 2 vary by a factor 16. The theoretical errors appear to 
be well under control, at least for this observable. 

So we can conclude that, in the future, the possibility to accurately determine 
the gluon distribution in the photon relies on the possibility to improve the experi- 
mental systematic errors. In the future also statistics will be larger and higher Et 
regions will be accessible. The corresponding data will put more constraints on the 
parton distributions. 
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